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ABSTRACT 

Continued studies were made of bluff-body separation and its 
effect upon the response of an aeroelastic model for a range of tunnel 
air speeds and model cross sections. The motivation for this work 
was the application to the antenna-element vibration problem as posed 
in the various Wullenwebber type antenna installations operated by the 
United States Navy throughout the world. 

A retrofit configuration change was devised to significantly 
decrease the wind-induced response of the long horizontal "boom boards". 
Although flow separation Mt ewiaccsieret the character of the aerodynamic 
forcing function has been modified such that the energy spectra have 
apparently been broadened from a single dominant peak into many weaker 
peaks. Thus the effective total response is reduced. 

The digital data processing has been extended in scope due to the 
repeatability in signal analysis. Estimates have been made of the "'two=— 
dimensional"' random input forcing function using the output response in 
conjunction with the aeroelastic-model transfer function. Reasonably 
good results were obtained and an estimate was made of the root-mean- 


square value of the unsteady fluctuating air load. 
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TABLE OF SYMBOLS 


D - model depth 

f - frequency (cycles per second) 

fT - natural frequency 

H - model frontal width 

i - mass moment of inertia of model relative to pivot axis 
ee - length of spring rod (inches) 

m - mass (slugs) 

q - test section dynamic pressure, 5 v2 (psf) 

RAY - ratio of glove radius to length of side of square 
Re - Reynolds number = Vh/y 

S - Strouhal number 

V - velocity 


Z(iu)) - complex impedance 


8 - damping ratio, fraction of critical 

: : 3 
e - density of test section air, (slugs/ft ) 
WW - kinematic viscosity 


Note: Symbols associated with a particular derivation are defined when 


first used in text. 





i: Introduction. 

A bluff body may be defined from an aerodynamic sense as a body whose 
drag forces are primarily due to pressure or form drag. The airflow does 
not progress smoothly around the body, but instead separates randomly and 
defines a turbulent flow region. As a consequence of the flow separation, 
fluctuating pressures or air loads are generated which can be resolved 
into unsteady air load components in both the lateral (lift) and stream- 
wise (drag) directions. In general, the fluctuating lift forces dominate 
and give rise to the most significant vibration modes. 

Physical structures are elastic and consequently will respond to 
dynamic air loads that are wind generated. Examples of the response of 
bluff bodies to unsteady air loads are frequently observed in every day 
events, examples of which would include the galloping transmission line 
problem, vibrating telephone wires, smokestack motion, ground-wind loads 
on missiles, and vibrating-antenna elements. There are many other 
situations of structural shapes responding to wind-induced air loads. 

It was the purpose of this study to investigate the bluff-body 
problem with particular application to an important element of a large 
antenna (Wullenwebberittype) system that is currently in use by the 
United States Navy. The basic arrangement includes a number of towers 
located in a circular pattern. Figure 1 shows a pair of towers supporting 
a horizontal boom board. The boom boards extend around the entire circle, 
and serve the function of supporting vertical reflecting wires. The boom 
boards are 28 to 56 feet long (depending on the installation) and are 
approximately 100 feet above the ground. There have been instances 
reported of large amplitude motion of boom boards which were directly 
exposed to the airflow, and even some instances of component failures. 

An important consideration in this study is that the cross sections 
used for the boom boards are not the same for every installation since it 
is the option of each contractor to perform the detail design. Figure 2 
is a sketch showing the basic cross sections used for the boom boards. 
Only one installation (Winter Harbor, Maine) employs the circular cross 
section. There are three basic types of rectangular cross sections. The 
solid laminated boom board is used in three general sizes. The hollow 
square and the inverted "U'"' are lighter per unit cross section area, but 


lack the strength of solid construction. The inverted "U" does have 
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short lengths of wood across the bottom which support the vertical wires 
and provide structural support to maintain the shape. 

In order to fulfill the objectives of this study, it was desirable 
CO. 

(1) Gain an insight into the factors influencing the response of a 

Simulated aeroelastic bluff body. 

(2) Determine possible fixes suitable for application to the actual 
antenna installation. 

(3) Attempt to estimate the character and magnitude of the aerodynamic 
input. 

This work has been sponsored by the Structures Division of the 
United States Naval Civil Engineering Laboratory (USNCEL), Port Hueneme, 
California under work request WR 600028, dated 23 August 1965. 

The investigation was conducted with Lieutenant R. H. Jesberg, 

U. S. N. Tests were conducted in the West Coast Research Corporation 
(W.R.C.) low-speed wind tunnel located in Building 234 at the United 
States Naval Postgraduate School, Monterey, California. For detailed 
information on equipment and data handling procedures used in this study, 
see "Data Handling Procedures in Support of the Study of Aeroelastic 
Problems in an Antenna Vibration Project" by Lieutenant R. H. Jesberg, 

DO oa etl dual. 

Sincere appreciation is extended to faculty and staff of the 
Department of Aeronautics for assistance and support which was vital to 
the success of this investigation. The author is indebted to Associate 


Professor L. V. Schmidt for his invaluable counseling and guidance. 
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Za Test Description. 

An aeroelastic bluff body model was installed to vertically span 
the 3.5 x 5.0-foot test section of the W.R.C. low-speed wind tunnel. The 
wind tunnel is of the continuous flow type, and has an atmospheric vented 
test section. Test parameters included tunnel dynamic pressure (q), 
model cross-section geometry, and model natural frequency (f,). 

The cross section of the model was altered by fitting various shapes 
of balsa wood on an aluminum I beam. The low density of balsa wood made 
it possible to make changes in cross-section geometry without signif- 
icantly altering the weight, center of gravity, and mass moment of inertia 
of the. complete model. 

The aluminum I beam was flexurally supported at the top (outside of 
the test section) and could move laterally in response to the unsteady- 
lift loads. A drag constraint was provided at the lower end by a simple 
drag cable support. The natural frequency of the lateral motion was 
controlled by adjusting the length of a spring rod (2,) attached at the 
base of the I beam. Figure 3 presents the variation of model resonant 
frequency with spring rod length. Extreme care was used in the design 
of model components in order to assure repeatability of results and the 
avoidance of model slippage. 

For information on the complete design of the aeroelastic model used 
to represent the antenna boom boards, see Jesberg [7]. Included in [7] is 
a description of the displacement transducer used to measure lateral 
motion at the model base, the data handling and tape recording systems, 
the method used to convert analog to digital data, and the digital com- 
puter program employed to convert voltage time history into auto-correla- 
tion and power spectral density forms. 

A listing of the run schedule is provided in Table I. Runs 22 
through 45 were part of the first phase of testing and took place during 
the period of 22 October to 2 November 1965. The objectives of the first 


phase may be stated as: 


RUN OBJECTIVE 
ZORto: 27 Familiarization and check out of test set-up with 


square and circular cylinders. 


28 to 30 Check the effect of rounding the corners on a square. 
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CHES clo Sis) Investigation of miscellaneous configurations. 

Se Omcu Check the effect of rounding the corners on a square. 

The natural frequencies used in phase one were 36, 42, and 82 cycles 
per second (cps) which corresponded to beam lengths of i. =O. OR 
and 3.0 inches respectively. Upon completion of the first phase of 
testing, the data that had been recorded on magnetic tape were reexamined. 
Frequencies of all indicial response checks and tunnel test data were 
confirmed by plotting a portion of the tape data playback onto an ink 
recording oscillograph. 

Following phase one, modifications were made in order to improve 
data acquisition and reliability. The spring rod support bracket was 
redesigned for ease of setting spring rod length and for improvement of 
spring rod constraint. The displacement sensor was relocated to the 
centerline of the aluminum I beam. In addition, the spring rod was case 
hardened. 

The second phase of these investigations began with run 46 on 24 
January and continued through run 84 on 25 February 1966. The objectives 


for this phase of the program were: 


RUN OBJECTIVE 
46 to 51 Check the effects of modifications and compare 


results with first phase data. 

Zeon 62 Check shroud for retrofit on currently installed 

boom boards. 

63 to 74 Check scaling laws and effects of H/D ratio. 

Jomo // Check feasibility of using a shroud on flanges 

of the circular boom Boel 

78 to 84 Investigate the forcing function. 

The natural frequencies used during this phase were different due 
to the equipment modifications that took place between the first and 
second phases. The three basic natural frequencies used were 55, 70, and 
89 cps, which corresponded to Re values of 6.5, 5.0, and 3.0 inches 


respectively. 
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Dic Scaling. 

Accurate scaling of a wind-tunnel model is a fundamental considera- 
tion in establishing an experimental program. Important variables in- 
volved in an aeroelastic simulation include fluid density, viscosity, 
velocity, model vibration frequency, stiffness, and mass distribution. 
Since the model was restricted to a single degree of freedom (in the 
lateral direction), the moment of inertia did not have to be scaled. 

All essential scaling variables are included in either Reynolds 
number (Re = VHA’) or Strouhal number (S = fH/V). The Reynolds number 
of the model and prototype were matched to assure that the aerodynamic- 
forcing functions corresponded. All models were operated at values of 
Reynolds number between 250,000 and 750,000. The full-scale antenna 
element with a frontal width of one foot would have a Reynolds number of 
700,000 for a wind velocity condition of 75 miles per hour (mph) at sea 
level. 

Strouhal number is the dimensionless vortex-shedding frequency of 
the air flow and a match of this number assures a correspondence between 
the model and the prototype dynamic response. In general the prototype- 
antenna elements have a natural frequency of about five cycles per second; 
and as suggested by USNCEL personnel, a dynamic matching should correspond 
to a wind velocity of about 75 mph on a beam with approximately one-foot 
frontal width. At the design point, the fundamental structural Strouhal 
number is about 0.05. Since the model indicated large amplitude responses 
at S = 0.14 and multiples above, it is unlikely that the boom-board 
natural frequencies could act to reinforce the amplitude response at wind 
velocities above about 30 mph. 

It is customary to present wind tunnel results in the form of dimen- 
Sionless coefficients, i.e., lift coefficient (C, ); drag coefficient 
(C)), etc. In contrast the dynamic responses reported here are not 
normalized and in general it will be necessary to make relative compari- 
sons to see improvements or changes due to configuration revisions. 

It is interesting to check for dynamic matching since a good indi- 
cation of dynamic model characteristics can be obtained. Figure 4 has 
been arranged to show the lateral response of the model for the case of a 
Square cross section. Parameters shown in Figure 4 include model size 


(H) and spring rod length (Q,). A casual glance shows that shortening 
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the spring rod reduces the dynamic response, but decreasing the frontal 
width does not consistently decrease the motion amplitude. Therefore, 
comparisons can be made only when taking into account the scaling 
parameters. 

It can be shown that the root-mean-square (R.M.S.) of the lateral 
response is proportional to the expression below for a single degree of 


freedom model: 


= = = si i (Ff) As om 
‘Ka a ie ane | ( (2rrf.)4 £1-(#F¥P +4 & ‘ears (1) 
where: 
y - lateral response at the lower end of the model 
H - beam frontal width 
L - beam exposed frontal length 
IL - system mass moment of inertia about the flexure pivot 
f - frequency, cps 
a mean square of the fluctuating lift, assumed to be uniform over 


model length 
d;- Spectral density of input aerodynamic forcing function 


Now by a basic definition, 


Cc = ig OAs = © ops As (2) 


Equation 2 provides a scaling relation when Dis expressed in terms of 
density with respect to either frequency or Strouhal number. Since 


f/f, ~ S/S at a given operating condition, the scaling law is: 


\e = J Gr, wal fox" (3) 
V¥ VCR ete \ fn 


This law holds providing that I, =I, and ‘il Di (sy J, a Pon ee 


If the Reynolds number effect can be neglected and if there were no 
aeroelastic feedback dependence, then one could assume a = J ea 
On this basis equation 3 will simplify further to: 


T se or Hi ‘depot (4) 


— SS 
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Comparison again will be made using the data on Figure 4 under the 
assumption of matching Strouhal numbers. The steady-state drag data as 
shown on Figure 8 leady to the assumption that Reynolds number effects 
are relatively unimportant for a sharp-edged square or rectangular cross 
section. The starting point for the comparison is from run 58 at q = 20 
psf. The reference values are listed in Table II, part A. In making a 
ratio comparison, the response with the smallest displacement value 
(case c) will be used as a reference. On this basis, calculations were 
made and are listed in Table II, part B. It will be noted that the 
ratios of R.M.S. displacements do not agree with the ratios based upon 
q, H, and f but in all cases the displacement ratios were larger. A 
glance at equations 3 and 4 immediately points out that the assumption of 
constant fe as assumed in equation 4 is probably not valid. Column 3 
in Table II, part B is the ratio showing the discrepancy and has the 
physical interpretation of being the increase in Cz due to aeroelastic 
feedback or dynamic motion. 

Figure 5 is a plot showing the apparent increase in fluctuating lift 
due to motion. The abscissa represents lateral motion normalized with 
respect to frontal width. The data points are the values obtained from 
Table II, part B and all values are in reference to an elastic condition. 
A simple rotation of the curve, so that with zero motion the ordinate 
has a value of unity, provides an estimate of the lift increase due to 
aeroelastic feedback. It should be noted that in estimating this curve, 
no information has been obtained as to the magnitude of the fluctuating 
Life. 

A later section of this paper will deal with estimating the 
fluctuating lift and an account will be given to aeroelastic feedback 
again. In the discussion to follow concerning cross sections applicable 


to the antenna problem, only relative comparisons will be made. 


LZ 


4. Antenna Applications. 
H/D Ratio 
All but one of the boom-board cross sections for the Wullenweber 
antennas are rectangular in shape. The H/D ratio varies from one to 


two and a half. 


Le—-D —> 
WIND x t 
= H 

I 


A series of runs were conducted in order to demonstrate the influ- 
ence of varying D while holding H constant. The three beams selected for 


this evaluation were: 


(Ll) 6.0 x 3.0-inch rectangle (H/D = 2.0) 
(2) 6.0 x 4.0-inch rectangle (H/D = 1.5) 
(3) 5.7 x 5.7-inch square (H/D = 1.0) 


The small error introduced when using the 5./-inch square was considered 
acceptable since the values used for comparison were corrected to similar 
values of Strouhal number. 

The responses for natural frequencies of 89 and 70 cycles per second 
are shown on Figures 6 and 7. Table III is the numerical results of the 
H/D ratio study. It will be noted that the amplitude response is approxi- 
mately inversely proportional to the H/D ratio (holding H constant). 

The physical meaning of these results is that the fluctuating 
pressures acting over the lateral faces of the bluff body are primarily 
controlled by the corners and not the width of the lateral faces. 

Corner Radii 

A frequent statement made for circular cylinders is that the unsteady 
lateral air load is of the same order of magnitude as the steady state 
drag while the unsteady drag is an order of magnitude less. Using this 
as a clue, insight into the behavior of circular, square, and rectangular 
cross sections may be gained by looking at Figure 8 which is a summary of 
a drag study made by Delany and Sorensen[3]. The square cylinder does 
not show a drag transition trait as for a circular cylinder for Reynolds 


6 
number values well beyond 1xl0 . However, placing a radius on the corner 
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edges does bring a transition trait into being. Since the circular-cross- 
sectioned aeroelastic model was relatively more docile than the square 
cross-section,,the material of Delany and Sorenson provided the motiva- 
ten tO try Corner radii. 

Figure 9 shows a series of cross sections that gradually transform 
the square cylinder to a circular cylinder. The corners of the 5.7-inch 
square were successively rounded in one-half-inch increments of radius up 
to and including a two-inch radius. This series of cross sections trans= 
formed the square cylinder into a circular cylinder in definable steps. 

It was anticipated that rounding the corners of a square would produce 
responses that fall between the responses of the square and circular 
cylinders. 

Typical results obtained are the responses for a model with a three 
inch spring rod length as shown in Figure 10. The curve farthest to the 
left, is the response for a square and it will be noted that the response 
rapidly exceeds the limit of the graph. The addition of one-half-inch 
radius corners (R = 0.5 inch) on the 5.7-inch square (H = 5.7 inch) 
results in some improvement (R/H = 0.0877). The one-inch radius 
(R/H 
(R/H 


0.176) shows marked improvement. The one and a half-inch radius 


0.263) exhibited a low amplitude response characteristic and the 
two-inch radius (R/H = 0.351) was comparable to the circular cylinder. 

The advantages that can be obtained by rounding the corners of the 
boom boards are obvious when looking at Figure 10. The solid laminated 
boom boards would be particularly well adapted to this fix. However the 
boom boards with hollow centers would be seriously weakened if the corners 
were rounded. 

R/H = 0.263 is suggested for the construction of new boom boards. 
However, a corner radius ratio of R/H = 0.176 would provide a considerable 
improvement over the square corners which are now being used. 

Gloves 

Rounding the corners of a rectangular beam provided well defined 
decreases in lateral amplitude response of the beam. There are two 
reasons that discourage the use of this fix on the currently installed 
boom boards. First, all boom boards that are hollow in the center would 
be seriously weakened if the corners were rounded. ‘<8econdly, the boom 


boards are physically located 100 feet above the ground. It would be 
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easier to install some additional piece of material to a beam than remove 
part of the beam. 

A shroud on the corners of a square was selected as a possible fix. 
Figure 11 is a sketch of this shroud or glove as it will be referred to 
in this paper. From the results of an earlier series of tunnel runs 
where the corners of a square were rounded in successive one-half-inch 
increments, a ratio R/H = 0.22 was selected for the glove. The square 
with a two-inch radius (R/H = 0.351) had the most effect in reducing the 
dynamic response when compared with the basic square cross section. 
However, a glove with an R/H ratio of 0.351 would be relatively large, 
thus greatly increasing the weight and frontal area of the beam. A 
compromise was made in selecting a glove radius so that undue penalties 
in frontal area and weight addition would not take place. whiles still 
providing the benefits of reducing dynamic response of the boom boards. 

Two sets of gloves (R = 0.875 inches and R = 1.25 inches) were 
constructed to be used with two squares (H = 4.0 inches and H = 5.7 
inches). This provided four possible configurations with gloves to check 


for possible scaling laws and to give some indication of a desirable R/H 


ratio. 
(1) 5.7-inch square with 1.25 +inth gloves (R/H = 0.22) 
(2) 5.7-inch square with 0.875-inch gloves (R/H = 0.15) 
(3) 4.0-inch square with 1.25 -inch gloves (R/H = 0.31) 
(4) 4.0-inch square with 0.875-inch gloves (R/H = 0.22) 


Configuration (2) would indicate if the R/H ratio selected were too large. 
Configuration (3) would likewise indicate if the R/H ratio were too small. 
All four configurations have been plotted in Figure 12 for a spring rod 
length of six and one half inches (ft, = 55 cps). All curves exhibited 
considerable improvement over the response of the square cross section. 
Figures 13 and 14 show the response of the 5.7-inch squares with and 
without gloves. A ratio of R/H = 0.22 has a clear advantage of a lower 
response below a dynamic pressure of 30 pounds per square foot. When a 
natural frequency of 55 cps was used the R/H = 0.15 demonstrated a larger 
peaked displacement at 10 pounds per square foot dynamic pressure. 
Figures 15 and 16 show the response of the four-inch square with and 
without gloves. R/H = 0.31 and R/H = 0.22 are relatively the same 


amplitude at f =/0 cps (0, = 5.0 inches). However, at fo = 55 cps 
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(RQ, = 6.5 inches) R/H = 0.22 has a definite advantage below a dynamic 
pressure of 30 pounds per square foot. 

It is clear that the best radius for rounding the corners of a 
square is definitely not the best radius to be used for gloves on that 
same square. The ratio R/H when using gloves should be about 0.22. 

Circular Cylinder 

A positive advantage was derived by shrouding the sharp corners of a 
square. Therefore a brief look was taken at shrouding spanwise flanges 
located at the maximum width points (@=2£96) on the circular cylinder. 

The recommended fix for the circular cylinder by Hite and Breckon 
[2] was the rotation of the beam so that the flanges would be located at 
the stream stagnation points. Utilization of this fix in the field was 
infeasible for various physical reasons. 

Hite and Breckon demonstrated that the flanges when located at G=t9O° 
were responsible for the large oscillations of the circular cylinder. 
Without the flanges the flow separated randomly aft of the 90 degree 
position (present location of the flanges on the circular boom board). 
With the flanges installed, separation of the flow was stabilized to occur 
at the flanges, concentrating the energy of the lateral force at a 
dominant Strouhal number. The concentration of random energy in the 
frequency spectra results in a violent response. Therefore if the 
flanges can not be removed or rotated to stream stagnation points, the 
sharp edge effect of the flange should be eliminated. 

Shrouding of the boom-board flange was expected to produce desir- 
able responses. A one-inch radius shroud was selected, which was a ratio 
of shroud radius to beam diameter of 0.198 (R/H = 0.198). The shroud was 
extended laterally (into the stream flow) only far enough to cover the 
flange. This was a distance of only one quarter of an inch on the model, 
thus the cylinder retained a relatively smooth shape as is illustrated in 
Figure 17. The shroud was glued to the front half of the beam to provide 
the smoothest possible surface for the air flow. 

The responses obtained are plotted on Figures 18 and 19. With a 
spring lever arm of three inches, the response of the possible fix 
exceeded the response of a square (approximately the result of a flange 
at 90 degrees) at a dynamic pressure of 15 pounds per square foot. The 


proposed fix will apparently also exceed the square's response at about 


Zl. 


10 pounds per square foot when a spring rod length of six and one-half 
inches is used. A larger shroud radius relative to beam diameter might 
give better results. Another solution might be to blend the shroud 
fairing smoothly into the cylinder contour without discontinuity in 
surface slope. 

Further expenditure of effort was not considered justified because 
the circular boom boards at Winter Harbor, Maine have suffered further 
deterioration. Replacement of the circular cylinder with rectangular- 


cross-section boom boards is under consideration. 


Ze 


Se Estimate of Unsteady Lift. 

The last six runs (78 to 84) of the program were conducted in order 
to gain some knowledge of the lateral forcing function. Both the co- 
efficient of lift and its spectral distribution were estimated. A sharp 
edged configuration was avoided for this evaluation so that a full set of 
data could be obtained without restriction on amplitude response. For 
this reason, the 5.7-inch square with one-inch corner radii (R/H = 0.176) 
was selected. The natural frequency of the beam was varied from 60 to 
89 cps by adjusting the spring rod length Oh) in one-half-inch incre- 
ments from six to three inches. Data were obtained at tunnel dynamic 
pressures up to 40 psf using intervals of two psf. 

The amplitude responses (recorded on magnetic tape) were digitized 
to permit the use of a high-speed digital computer [7 ]. Utilizing the 
auto-correlation function, the power spectral density (P.S.D.) of the 
system output response was computed by the Tukey-Blackman method [1]. 
Sample power spectral densities of the output are shown on Figures 20 
and 21. These curves are computer plotted and represent the complete 
set of £. values at dynamic pressures of eight and 40 psf respectively. 
It will be noted on Figure 20 that there are energy peaks at 24, 50, and 
80 cps. These correspond to Strouhal numbers of 0.14, 0.28, and 0.45 
respectively. The large energy peak at S = 0.14 was evident for values of 
q up to about 15 psf (Re = 340,000). The characteristics shown on 
Figure 21 at q = 40 psf (Re = 560,000) indicate energy peaks at f = 41, 
80, and 123 cps which correspond to Strouhal numbers of 0.11, 0.21, and 
0.32 respectively, and in general the majority of energy is at the higher 
Strouhal numbers. These traits were shared by the major portion of out- 
put response data when Reynolds number was above 340,000. 

Computation of the power spectral density (P.S.D.) of the output 
from recorded data opens the door to investigation of the forcing 


fugetion. Dots) the P.S.D. of the output response was normalized. 


f Po (4) df =1.0 (5) 


The upper limit of i25 cycles per second assumes that there is no signi- 
ficant contribution to Ds (4) above that frequency. The accuracy of 
this assumption was checked and found to be correct for the calculations 
used in this report [7]. 


3 (5) is the P.S.D. of the input forcing function, (a distributed 
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force acting along the length of the model, but which may be treated as 
a concentrated load acting at the mid-span of the model). Di. (+) is 


related to DoF) through the following relation. 


e 
Oi = |ZzaHl ow (6) 
& (it) is the transfer function which is dependent only on the model 


being tested. 


lzasl’={Li.0-(+) J +4 ¢ (£)} (en fn)” (7) 

f is the frequency in cycles per second and rT is the natural 
frequency of the beam. B is the damping ratio of the beam. It was 
calculated by log decrement procedures on indicial-response checks. Due 
to the difficulty of reading the amplitudes of successive peaks with 
sufficient accuracy, an average B of 0.018 was used for computation of 
Di (4) . This value was obtained from eleven log decrement calculations. 
The damping ratio was sensitive to spring rod clamping. Various values 
of 5) (0.018 to 0.027) were used in calculating QD; (+) to check the 
SLicemeO mE Aw 0apen centeincrease in @ resulted in less than one 
per cent change in Q,; (4) 

The forcing function may be found by considering the geometry of 


the test set up in the wind tunnel as sketched below: 


Values used were: 








See Crive 
SIMPLE 
SUPPORT 
PIVOT 


Xq = 4.0 inches 


X,, = 46.0 inches 
MODEL 


= 257 Oeinehes 
X4y = 54.0 inches 
I = 35.73. slug-in- 


DISPLACEMENT — 

PICKUP ae sl 

rn 

—spRING CONSTRAINT 


The equation of motion for harmonic excitation may be expressed as: 


T6+céeé+kxe = F LE Join wt (8) 
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Replacing the angular coordinate by the displacement at the transducer 


yields: 
¥ = X40 
E Xe Apt Xa— 
SY + Kee y = ‘a =o | sin wt 
v, : F Xa [ An+X 
Yen rae Con +unY =F [ Ae an et 
A FE Xa Xp + Xa 
let M = —+ L 2 ] 
Mm 
Nz IZ cia] 
as — 2 
my =F lZzciwl 
If the mean square of the displacement is ve , then the mean square of 


the applied load is: 


= xa Xp + Ka 2. 
ha =- =| = = (M.5.)( ¥4) 


pies. 
where (M.S.) = i, De (+) | Z is) I" AS 


° 


Ae 


C= De (epee 


and ay _ T= 
\ Ch a % H (Xo-Xa) 
a M.S.) CY# ah (2L)_ 
therefore (Gow ; H Xa ( Ke “X29 (9) 


‘ioe is the root-mean-square of the coefficient of lift. 

A high-speed digital computer was utilized to calculate the root- 
mean-square coefficient of lift for each of the seven natural frequen- 
cies of the beam and at each dynamic pressure used. Figure 22 is the 
plot of the individual Ara and Figure 23 is a plot of the average Ice 
Both plots clearly indicate the occurrence of a transition at 14-pounds 
per square foot dynamic pressure (approximate Reynolds number of 
Ses a 107) where the flow exhibits a change to supercritical flow. 

All curves are essentially parallel and with only minor exceptions the 


beams with lower natur3l frequencies have larger coeffictents of lift. 


Z 


The consistent pattern of these plots suggests that aerodynamic forces 
may be induced by the motion of the models through an aeroelastic-feed- 
back influence. 

The P.S.D. of the input was also machine plotted using the digital 
computer to handle the transfer function type calculations. Figures 24 
and 25 represent estimates for subcritical and supercritical flow and 
show good agreement for all values of model resonant frequency With 
respect to the occurrence of energy peaks. 

The peaks of the power spectral density curves of the exciting 
force occurred at a constant Strouhal number. For dynamic pressure up 
to and including 12 pounds per square foot (Re = ees the vortex 
shedding occurred at a Strouhal frequency of 0.140. At and above 16 
pounds per square foot (Re = jG the Strouhal frequency was 0.32. 
At 14 pounds per square foot (Re = aa two peaks were present, 
one at a Strouhal number of 0.14 and the other at 0.32. The beams with 
natural frequencies below 75 cycles perssecond (70 65, and 60 cycles 
per second) responded at a Strouhal number of 0.32. 

The 14 pounds per square foot reading was obviously in the flow 
transition region. Supercritical flow was experienced at a dynamic 
pressure of 16 pounds per square foot (Re = 4S ell oy Dr. Fung found 
that the critical Reynolds number for his circular cylinder was 
Re = 3.0x10° [6]. 

The aeroelastic-feedback term was further investigated by consid- 
ering ie results of Figure 22 at dynamic pressures of 10 and 30 psf 
which correspond to Reynolds numbers of 280,000 and 480,000 respectively. 
The fluctuating lift values are cross plotted on Figure 26 as a function 
of the normalized lateral response. It will be noted that both flow condi- 
tions have approximately linear dependence upon motion, and the sub- 
critical flow appears to extrapolate to a value of Vc = 0.35 for a 
rigid model while the corresponding supercritical flow value is about 
0.075. A graphical procedure has been used on Figure 26 to obtain an 
aeroelastic gain correction for the subcritical flow case, normalized 
such that the gain value is zero for the rigid case. The large influence 
of motion on fluctuating loads is quite surprising, especially when one 
notes that lateral displacementy(R, 5.) on the order of.one hair threk- 


ness more than doubles the value of fluctuating load. 
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A simple model may be conjectured to describe the build-up of an 
aeroelastic feedback term. The model assumes that fluctuating lift is 
developed due to the instantaneous response of a rotating drag vector, 
with the rotation being caused by the motion-induced angle of attack 
changes. Consider a bluff body moving laterally in rectilinear flow as 


shown in the sketch below: 


If y = o sinwt, 
then del. =Y= Ye sin wt 
hence Unax = ert Ye 


—— i= 
and U" = eyapee cr y* 


Therefore the mean-induced angle is given by: 





oa ee ee te 
a yo le ae 
But if there exists a dominant Strouhal number that defines the frequency, 


then f = VS/H. 


—— = 
Hence 4a eile o* (Zz ) 


Lf we assume that Cr (+) = (+) Cp STEADY 


then Vcr =Ve'rs ve CO ay (10) 


and “A fay = VE WS Costcacy i) 


For the square cylinder with rounded corners under the condition of 


subcritical flow, we may assume S = 0.14 and Ch = 1.0. Hence an estimate 


Za 


of the aeroelastic effect would be 


ary en EY =Ve TY (0.14) (1.0) = 6.2 %* lo! 


The actual slope estimated from Figure 26 is 5 Salen a value 
approximately four orders of magnitude greater than that predicted by 
the simple model. The large discrepancy would indicate that either the 
experimental procedures and estimates presented here are suspect of that 
the simple model is completely inadequate. 

The existence of fluctuating loads on a rigid model seems reasonable 
Since the unsteady pressure exist in separated flow, independent of motion. 
The dependence of these pressures and their instantaneous integrated 


value (net fluctuating force) upon motion is not clear, and further study 


certainly seems worthwhile. 
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6. Summary. 

The large responses of the square-cross-section boom boards were 
sharply reduced by two different modifications to the beams. Rounding 
the corners on solid boom boards at the time of construction appeared to 
be a desirable solution. A ratio of corner radius to frontal width of 
6.263 provides a response that is near the level of a circular cylinder 
for all wind velocities (up to 171 miles per hour) that may be encountered. 

A retrofit modification was tested for use on boom boards that are 
presently installed as part of the Wullenwebber Antenna system. This fix 
consisted of shaped shrouds which would be attached over the corners of 
the boom boards. The desirable ratio of shroud radius to the length of 
the side of the square should be Abomt(0.22ic C.01.. 

A shroud for the flange on the circular cylinder was tested, but no 
appreciable improvement was noted. A different size and or shape of 
shroud might produce a marked reduction in amplitude response, but fur- 
ther studies in this direction were not made. 

Tunnel tests demonstrated a correlation between models of different 
H/D ratios (ratio of frontal width to depth of model) in the Reynolds 
number supercritical flow region. The first order components of the 
amplitude responses were found to be inversely proportional to the H/D 
ratios (when holding frontal width constant), thereby showing that the 
fluctuating pressures developed on the lateral faces are not too sensi- 
tive to the length of that face. The dominant factor appears to be 
corner radius ratio. 

The use of precision instrumentation, magnetic recording of data, 
and data processing by analog and digital computers permitted investi- 
gation of the forcing function's power spectral density and lateral 
root-mean-square coefficient of lift. The coefficients of lift were 
estimated for the square with rounded corners and were compared with 
previously published information on circular cylinders. The dependence 
of unsteady lift upon motion strongly suggests the influence of an 
aeroelastic feedback term. Further investigations should be conducted 
in this area with different size models to confirm and/or isolate the 


aeroelastic feedback component. 
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IMaAbLe T 


RUN SCHEDULE 


Date Run Gonfiguration q(PSF ) R Data 
10-22 20  Gire Cyl 20 8.88 F,A,D,T 
10-26 21 Circ Cyl 5-10 8.88 F,A,D 
22. Cire Cyl 5-40 8.88 F,A,D,T 
Ze Cire Cyl 5-40 6 .00 
24 Gulige (Gy 5-40 $) 1G 
10-27 25 Sqr Cyl 5-20 200 
| 26 Sqr Cyl 5-12 5.00 
27 Sqr Cyl 3- 8 6.50 
10-29 28 5.7 Sqr & .50 RC 3-10 6.50 
29 5-15 5.00 
30 5-30 3.00 
om Front-5.7 Sqr & .50 RC 5-25 3706 
Rear-Rect Cyl 
10-30 32 u 3-10 6.50 
35 Rect Cyl 3- 8 6.50 
34 | 5-15 5.00 
35 5-20 8200 
36 Front-Rect Cyl 5-15 300 
Rear-5./ Sqr & .50 RC 
11-1 37 Sqr Cyl & 1.00 RC 5-40 3.00 
38 | 5.00 
39 6.50 
40 5.7 Sqr & 1.50 RC 6.50 
41 5-40 5.00 
42 3.00 
pl - 2a 43 Do esdr -& 2200) RE 3.00 
| 44 6250 
45 STAG, 
1-24 46 5.7 Sar & 2.00 RC 5=40 5.00 D 
47 Seyeogqr & IeZ5 Glove 3.00 D 
42 | 5.00 D 
49 6.50 D 


GE 


Date Run 


2-11 50 
sil 
2-14 52 
53 
54 
55 
56 
57 
2-15 58 
59 
60 
61 
62 
2-18 63 


Zoe 78 


Y WabYe ‘I 


RUN SCHEDULE (CONT'D. ) 


Configuration 


a eee 


aye 


Sf Fe. uu GS 


7 Sqr & 1.25 Glove 


-/ Sqr & 1.25 Glove 


-/ Sqr & .875 Glove 


| 


a oer e he 

i aod 

.0O Sqr 

.OO Sqr 

-0O Sgr 

nC SQray 675 


| 


-OO Sgr 1.25 Glove 


.00 x 4.00 Rect 


7 


.00 x 3.00 Rect 


Circ Cyl 1.00 R over 
flange 


a) 


: 


7 Sqr & 1.00 RC 


Y 
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q(PSF) 2 Data 
5 ,40 6.50 A,D 
3.00 A,D 
S700 FADE 
5.00 
6.50 
6750 
5700 
3.00 
5-20 500 
5-15 a. OC 
5-12 UG 
3- 8 6.50 
5-20 300 
5-40 o OOwee! .A,D, 1 
510, 
6.50 
> OO 
3.00 
6.50 
3- 8 6.50 
3-15 Sl: 
Siro) 3.00 
5eZ5 300 
5-20 3) 31010, 
Sz 6.50 
5-40 6.50 
5-60 500 
5<35 3.00 
5) TA0, D00°"° uP DAT 
Y 3.50 + 


fable: Ll. 


RUN SCHEDULE (CONT'D. ) 


Date Run’ Configuration q(PSF ) Q. Data 
2-25 80 5.7 Sqr & 1.00 RC 5-40 4.00 F,D,T 
81 4.50 
82 5.00 
83 D250 
84 6.00 

RC - Radius of Corners 


Cire - Circular 
Cyl - Cylinder 
Sqr - Square 

Rect - Rectangular 


Q,, - Length of Spring Rod 


Data Code: 
F - Frequency Response Measured (Tunnel-off) 
- Acceleration Readings Recorded 


- Displacement Readings Recorded 


3 Oo } 


- Data Recorded on Magnetic Tape 


All lengths in inches 


BS) 


Table II 


Scaling Study Calculations 
Square Cross Section 


Reference Reference Data 





B. Comparison 


(d):(c) 


be) Gee) 
(a):(c) 
(b):(c) 





34 


Table III 
H/D Calculations 


a) 
~ 
7) 


(10 


BEE 


oO WwW OF oO UV oY} CO UW O10 mM OF Oo mM OF oOo OM oO 


pee 


bee 


Ze 
1 
1 
2 
1 
1 
2 
1 
1 





r= KN 


[= 





* Amplitude ratio is defined relative to the amplitude of the bluff 
body with H/D = 1.50. 
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FIGURE {| 
ANTENNA INSTALLATION 
GENERAL LAYOUT 
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FIGURE e2 
TYPICAL BOOM BOARD 


CROSS SECTIONS 
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FIGURE o3 
SPRING BAR CALIBRATION 
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FIGURE 4 
DYNAMIC RESPONSE, SQUARE CYLINDER 
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DYNAMIC PRESSURE ,g (PSF) 


FOU 
SFreOGTS OF MOTION ON’ CiF fT, 
SQUARE CYLINDER 


p 


Data From 
| RUN to RUN 
©O 62 62 
A 60 " RUN 62 Is Elastic Ref. Value. 
O 61 ~ see TABLEI for Computation. 
V 58 : 
© 59 





oe Estimated Lift Increase 
ft Gain Rel. to Rigld Value 
Eactor I 


. Wa 


Lift Increase Relative to 


a Elastic Ref. Value (0) 
fi 
— O 


eo) 0.4 0.8 1.2 1.6 2.0 


4) y® “Normalized Lateral Motion 
H -RMS-( x 10°3) 
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FIGURE 6 
EFFECTS OF H/D RATIO, 2.23 INCHES 


RECTANGULAR CROSS SECTN.CYLS., Q=3.0in. 
SYM = RUN 


-p i. 
A 1.0 58 
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FIGURE 7 
EFFECTS OF H/D RATIO, 2,=5 INCHES 


RECTANGULAR CROSS-SECT. CYLS., &,=5.Oin. 
SYM = RUN 
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FIGURE 9 


MODEL CROSS SECTIONS 
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FIGURE I] 
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Figure 20 
P.S.D. OF OUTPUT AT q=8 PSF 


Notes: 


Square cylinder and corner 


radii, R/H=0.176 


Re= 2.50% 10° 


7 curves shown summarize {, 
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change from 6.0" to 3.0" 
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Figure 21 
P.S.D. OF OUTPUT AT q=40 PSF 
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Notes; | 
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Figure 24 
P.S.D. OF INPUT AT q=8 PSF 


3. 7 curves shown summarize g., 


change from 6.0" to 3.0" 


Notes; ; 
1. Square cylinder and corner 
radii, R/H= 0.176 
2. Re= 5.6 x 10° 
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Figure 25 
P.S.D. OF INPUT AT q=40 PSF 








Notes: 
Square cylinder and corner 
radii, R/H=0.176 


Re=2.5 x 10° 
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FIGURE 26 | 
AEROELASTIC EFFECTS ON UNSTEADY LIFT 


SQUARE CYL. +CORNER RADII; © =0.176 
SYM Re 
6 
© 0.28x10", 10 PSF 
A 0.48x10°, 30 « 
Note: RMS of C, Obtained 


from Digital Computer Estimates. TY = 
- 4 
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air load. 
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